Multiple scattering of waves by a pair of gravitationally stratified flux tubes 
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ABSTRACT 

We study the near-field coupling of a pair of flux tubes embedded in a gravitationally stratified 
environment. The mutual induction of the near- field jackets of the two flux tubes can considerably 
alter the scattering properties of the system, resulting in sizable changes in the magnitudes of 
scattering coefHcients and bizarre trends in the phases. The dominant length scale governing 
the induction zone turns out to be approximately half the horizontal wave length of the incident 
mode, a result that fits in quite pleasantly with extant theories of scattering. Higher-/? flux tubes 
are more strongly coupled than weaker ones, a consequence of the greater role that the near- 
field jacket modes play in the such tubes. We also comment on the importance of incorporating 
the effects of multiple scattering when studying the effects of mode absorption in plage and 
interpreting related scattering measurements. That the near-field plays such an important role 
in the scattering process lends encouragement to the eventual goal of observationally resolving 
sub- wavelength features of fiux tubes using techniques of helioseismology. 

Subject headings: Sun: helioseismology — Sun: interior — Sun: oscillations — waves — hydrodynamics 
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1. INTRODUCTION 

An outstanding issue in solar physics concerns the accurate constraining of the internal constitution of 
sunspots. Since we are unable directly image the interior, we study the solar acoustic wave field in and around 
sunspots and attempt to comprehend these observations through theories of wave inter actions. An exampl e 



of such an effort is the first putative detection of downflows underneath sunspo ts by Duvall et all (jl996l ). 



who analyzed the solar wave field using methods of time-distance helioseismology ( Duvall et al. 1993h . In the 



Sun, observations have almost always been more plentiful than theory. Of late, the importance of developing 
theoretical and computational methods to aid the interpretation of observations of solar magnetism has 
risen to the fore. A considerable body of computational wor k has recently focused on understanding the 



nature of wave intera ctions in magnetized environments (e.g. iKhomenko fc CoUados 



2006 



Cameron et al 



2007 : Hanasog^ 2008 ). Complementary to such efforts, we attempt here to develop further the theory of 
fiux tube related multiple scattering, probably ubiquitously present on the Sun but has, for the most part, 
been studiously ignored in the past due to the many challenges involved in modeling such interactions. On 
another front, it is important to place bounds on the degree of wave absorption and scattering by the plage 
for it tells us how much energy is transmitted to the corona and could help explain the co mplex frequency 
dependence of acoustic mode linewidths (e.g. Bogdan et al. 19961 : Hindman fc Jain 2008h . Because plage 
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comprises ensembles of compactly packed thin flux tubes, the interaction with waves possibly lies in the 
multiple scattering regime. 

In linear theory, when a wave encounters an anomaly of some sort, it is scattered at constant frequency, 
with the resultant wave field being broadly classified into near- and far-field components. The far-field 
consists of propagating modes that transport some fraction of the incident mode energy away from the 
scatterer. The near-field is more complex, comprising a number of non-propagating horizontally evanescent 
waves, that arise when the displacements of the anomaly due to the external wave buffeting cannot be 
matched by the set of eigenfunctions of the propagating modes. In the case of stratified fiux tubes in the 
Sun, the set of p-mode eigenfunctions at constant frequency is an incomplete basis, requiring a supplementary 
set, an uncountably infinite continuum in fact, of these evanescent near-field functions to complete the basis. 
The problem is exacerbated when the displacement eigenfunctions of the scatterer gain complexity. 

The mathematics required to address the near-field in the case of thin fiux tubes was set down by 
Boedan fc Calivl (|l995l) . who termed this non-propagating sheath of waves as an "acoustic jacket" that 
envelopes flux concentrations. The uncountable continuum of jacket modes arises as a consequence of an 
infinitely deep lower boundary, required in order to allow tube modes to disappear into the solar interior. 
Unfortunately, numerically computing the near-field jacket in this scenario is all but impossible because of 
various formidable inte grals in the equations. In order to arrive at an analogous but more tractable problem, 
Barnes fc Callvl pOOOT ) introduced an artifici al lower boundary, lead ing to a discrete and c ountably infinite 
number of near-field modes. More recently, iHanasoge et al.l (j2008l ) adopted the model of iBarnes fc Callv 



pood ) and estimated the magnitude of the near fiel d jacket and the sin gle scattering by an isolated thin flux 
tube. They attempted to model the observations of lDuvall et al.l (j2006[ ). who characterized the scattering of 
/ modes by magnetic flux elements. However, magnetic elements are known to consist of a number of tightly 
packed flux tubes, all likely within the near-fields of each other. Thus the single scattering assumption may 
not be entirely accurate when interpreting these measurements. 

The presence of a large body of theory to draw upon makes it easier to proceed towards an understanding 
of the importance of multiple scattering. When a pair of flux tubes lie in the proximity of each other, their 
near-fields communicate and depending on the separation, can dramatically alter the nature of the scatter. 
Thus accounting fo r mutual induct i on of near-fields is rather important when studying plage or other closely 
spaced scatterers. iBogdan &: Foxl ()l99l[ ). when considering a pair of fiux tubes at a series of separations 
in an unstratified medium, found evidence for three different scattering regimes that they termed multiple, 
coherent, and incoherent. The nomenclature points to differences in the degree of coupling between the two 
flux tubes, with the incoherent regi me no different from is olated body scattering and the multiple regime, 
substantially different. Subsequentlv iKeppens et al.l (|l994 ) studied ensembles of flux tubes and found that 
the degree of absorption was greater in "spaghetti" models than monoliths, pointing to a way of discerning 
the differences between the two. 

These efforts have been restricted to unstratified media, mainly due to the considerable mathematical 
complexity that stratification injects. With the addition of gravity, the external driver and the fiux tube 
displacement eigenfunctions assume distinct and more complicated forms, possibly destroying a number of 
resonances hitherto possible in the unstratified case. Furthermore, purely analytical techniques cease to be 
of utility, even when consideri ng single scattering, let alone its multiple counterpart. When studying mode 
mixing due to thin flux tubes, Hanasoge et al. ( 2008h were forced to apply a number of methods of linear 
algebra in order to reliably estimate scattering coefficients and the near-field jacket. On the other hand, 
multiple scattering is a bit more of a challenge since one must simultaneously determine the wave fields of 
a number of disparate scatterers, all the while keeping in mind that each scattering coefficient is defined 
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according to a coordinate system centered on that specific scatterer. Fortunately, the theoretical machinery 
to s tudy these sorts of prob l ems has been develop ed decades ago, in the context of fluid mechanics, by 
e.g., Kagemoto fc Yuel |l986l ); lLinton fc Evans! (fl99a' ). We adopt these methods in our calculations in order 
to determine the degree of mode mixing and scattering from a two-tube system. Please note that in the 
discussions below, the terms "near field" , "acoustic jacket" and "envelope of evanescent modes" are used 
interchangeably. The plan of this paper is a s follows. We describe th e stratification and basic aspects of the 
flux tube tube model in g2 The method of iKagemoto fc Yuj (|l986l ) in the context of interacting thin flux 
tubes is discussed in ij3l The scattering coefficients derived through the application of these techniques for 
the two-tube system for different incident modes are presented in 21 Finally we summarize and conclude in 



MODEL 



The background structure in this calculation, adapted from lBogdan et al.l (jl996f ): lHanasoge et al.l (j2008l ). 



is an adiabatically stratified, truncated polytrope with index m = 1.5, gravity g = —2.775 x 10'* cm s '^z, 



reference pressure pa — 1.21 x 10^ g cm" 
the pressure and density variations are given by. 



'"'^ s ^, and reference density /9o = 2.78xl0 ^g cm such that 



Piz) = Po 



zo 



m+l 



(1) 



and 



piz) = Po 



Zo 



(2) 



We utilize a right-handed cylindrical co-ordinate system in our calculations, with coordinates x — (r, 9, z) 
and corresponding unit vectors (f, 9, z). The photospheric leve l of the background mo del is at z = 0, with the 
upper boundary placed at a depth of zq = 392 km. Following Barnes &: Callv ( 2000[ ). we introduce a lower 
boundary at a depth of 98 Mm. The displacement potential ^'(x, t) describing the oscillation modes {t is 
time) is required to enforce zero Lagrangian pressure perturbation boundary conditions at both boundaries. 
This upper boundary condition is reflective in nature and therefore, po ssibly not very rea listic. The incoming 
p„-mode, a plane wave, which expanded in cyhndrical coordinates (e.g. lGizon et al has a displacement 

eigenfunction, ^E'inc, of the form: 



niO—ut] 



(3) 



where, 



The various symbols in equations ([3]) and ^ are: 
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uj the angular frequency of oscillation, s = — z/zi^, Jm(w), the Bessel fu nction of order m and argument 
w and MK^fiiw), the Whittaker function (e.g. IWhittaker fc Watson Il963h with indices k, n and argument 
w. The eigenvalue > and constant characterizing the mode are obtained through the procedure 
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described in appendix A of Hanasoge et al. ( 2008[) . The n = mode corresponds to the surface gravity or / 
mode, while n > represents the acoustic p„ mode. Note that the lower boundary results in a finite sized 
box and hence places a restriction on the number of p modes that can fit in this domain. The term Nn is 
the normalization constant for the mode, defined as 



9 



ds 



-1/2 



(6) 



The near-field eigenfunctions are also solutions to the same differential equation that governs the prop- 
agating modes: 



Qn{l^n \ s) — S /- ^ 



(7) 



As can be seen, the only difference between the form of the propagating and evanescent mode eigenfunctions 
is the fact that the roots are now imaginary. In order to determine the roots, we perform a high resolution 
search for eigenvalues k^'^ and constants C^'"^; the task is relatively easy for the propagating mode parameters 
but is reasonably difficult for the jacket modes because the eigenvalues may be very finely spaced (appendix 
A of .Hanasog e et al, 2008) . Subsequently, tables of the propagating and jacket modes are pre-computed at a 
range of frequencies. Computations of Whittaker functions over large parameter spaces are fairly non-trivial; 
we use a number of CERNLIB routines to accomplish all these tasks. 



2.1. FLUX TUBE 



Applying the approximations listed in §2 of iBogdan et al.l (|f 996l ). a thin flux tube carrying a magnetic 
flux of $/ = 3.88 X lO^^Mx, with constant plasma-/? everywhere inside the tube is embedded in the polytrope. 
The thin flux tube approximation. 



/87rp(s) 
1 + 13 ' 



7ri?2(s) 



(8) 



where b{s) and R{s) are the magnetic field and the radius of the tube at depth s, is shown to be accurat e 



1996h . 



to better than a percent in the truncated polytrope situated below z = —zq or s = 1 (jBogdan et al. 
Note that the magnetic flux associated with the tube is held constant - different values of /? therefore result 
in different b{s) and i?(s). The constant-/? property of the tube follows from the assumption that thin flux 
tubes are for all practical purposes in thermal and radiative equilibrium with the external medium. 



2.2. Oscillations of the tube 



For the cases addressed here, we treat only horizon tal kink motions of the flux tube {^±{uj,s)), caused 
by impinging m ~ ±1 modes (e.g. iBogdan et al.lll996( ). The m = ±1 modes affect the tube oscillations 
according to the differential equation: 



to zo + 



2gs 



9 



d 
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1 + 2(3 



/S) 2 
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dx 



(9) 



{l + 2l3){m + l)ds^ l + 2l3ds_ 

where x — rcosO. The scattered wave field is computed by matching the horizontal components of the 
motion of the flux tube to the external oscillation velocities. The manner in which this is accomplished is 
detailed in the following section. 
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METHOD 



Consider a system of randomly distributed flux tubes. The scattered wave field around tube i with the 
origin of the coordinate system located at the center of the upper boundary of the tube is given by: 



E 



np 



N 

E 

n—n p -f 1 



(10) 



where Cn are near-field eigenfunctions, $„ describe the propagating p-modes, Hm (x) and Km{x) are Hankel 
and K-Bessel functions of order m acting on argument x respectively, n p denotes the number of propagating 
mode eigenfunctions (a finite number due to the presence of the lower boundary) and the rest corresponding 
to the evanescent jacket modes. Also note that the m su mmation is truncated, sin ce thin flux tube theory 
applies only to interactions with |m| < 1 waves. Following Kagemoto fc Yue ( 19861 ). we write this in matrix 
form: 



where Ain is a vector (of size 3 x 1) of scattering coefficients for tube i and mode n. The matrix (size 
3 X , with being the number of points in the z grid) contains the partial wave expansions in terms of 
In particular, the elements of A. are: 



A ■ 

■^171 



a, 



On 



{n < np), 
{n > np), 



(12) 
(13) 
(14) 



where is the dth point along the s axis, and the indices c,d run from [1,3] and [Ij/i^], respectively. 
As described in the introductory section, the challenge in computing the wave fleld interactions lies in 
simultaneously solving for the scattering coefficients of all the tubes. Tubes that are placed sufficiently far 
from each other can only interact via the far-field propagating modes since the evanescent jacket has a spatial 
decay scale of a wave length or so. Thus the coupling is reasonably weak since the amplitude of the expanding 
far-field modal ring falls rapidly with distance from the scatterer. Stronger interactions occur when tubes lie 
within each other's near-field induction zones, for the evanescent modes of one influence the tube oscillations 
on the other and vice-versa. Fairly significant changes in the scattering cross-sections are a consequence 
of this phenomenon. Thus, to capture this effect, we must compute the wave field around each tube and 
project its influence on the neighboring flux elements. Since the wave field around each tube is written 
in a co-ordinate system centered along its axis, suitable co-ordinate transformations must be performed. 
For this purpose, we em ploy Graf's addition formulae , which transform cylindrical wave functions between 
co-ordinate systems fe.g. lAbramowitz fc Stegunlll964l ): 



H^Ji\Kn)e 



mi(0i--iii) _ 



E 



H'^^Xa{KRu)UKn)e'"^^- 



l+lil) 



(15) 



if™(fc,{rOe*'"(''-^") = K„,+d{K{Ru)Id{k'^n)e''^^-''+''^^\ (16) 

d — — CG 

where Im{x) is a Bessel-I function of order m acting on argument x. These equations show how to expand 
an outgoing scattered wave function from tube i in terms of the incident mode waves of tube I. Here Rn is 
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Transformation parameters 





,1 















X 



Fig. 1. — The transformation parameters of equations (fT5|) and (fT6|) shown graphically. The tubes i^l are 
separated by a distance Ru, with the line joining their centers inclined at angle to the x-axis of tube i. 
Note that ^u=tt - ^u. 

the distance between the centers and ju the angle that the line between the centers of the tubes subtends 
at the X-axis of the tube i (shown graphically in Figure [1]). This relation shows how the scattered wave field 
from i acts as an incident wave on I. Written in matrix form: 

= T^i^L (17) 

where T is the transformation operator that relates the scattered wave field of tube i to the resultant 
incident field on I. We use equations (fT5|) and (|16p to build T; the precise distribution of elements is listed 
in appendix 1X1 Note that the index n that appears in equation (|17p denotes a specific modal order, while 
the superscripts /, S represent the incident and scattered wave expansions. Recalling that equation (jlip 
contains the contributions of various near- and far-field waves and collecting all the n's, we have: 



Equation (fTS]) tells us how much the scattered wave field of tube i acts as an incident field on /. Summing 
up the contributions from the scattered wave fields of all the other flux tubes (except itself, of course) and 
the zeroth order incident wave, the total incident wave field at I is given by: 



TV 



= EUol'"+ E ^iT^u'^L] (19) 

- EK"+ E ^»Tr,UL (20) 

and (/)o|iri = o.fn'^hi^ where is a vector of coefficients representing the zeroth order incident wave mode of 
order n, acting on tube I. The determination of the amplitudes of the off-axis incident modes is described 
in appendix |B] Note that ain and Ain have the same matrix sizes; the former refers to the zeroth order 
incident wave whereas the latter contains th e net scattering terms. F inally, to close the equations, we use the 
diffraction transfer matrix approach of e.g. iKagemoto &: (|1986I ) to relate the total incident coefficients 
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to the scattered ones. In effect we generate a matrix B; that acts on the incident wave field at / and produces 
the scattering coefiicients. We know the total incident wave field at I: it is given by equation ^U\i . And 
since Ai is the vector of scattering coefficients as seen by the co-ordinate system centered on I, the following 
relation must hold: 



A, =B 



or, 



N 



(21) 



(22) 



where we have the following relations: 
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(23) 
(24) 

(25) 

(26) 

The matrix B/ is constructed as follows. The scattering into all other modes is computed for each incident 
mode {m,n). This includes both propagating p„ modes and near-field type incident waves (the Bessel-I 
functions). Thus B/ contains a full description of the scattering of any given incident mode (to, n) into any 
{m',n') for tube I. More details are described in appendix \X[ 

Finally, we discuss the means applied to study the simple case of two interacting thin flux tubes 1,2. 
Proceeding from equation (I22|) we have: 



Ai = Bi (ai +T^iA2) , 

A2 = B2 (aa + TfaAi) , 

Ai = Bi (ai+T^iB2 [a2- 

[I - BiT^iB2Tf2] Ai = Biai + BiTfiB2a2. 



TliAi]), 



(27) 
(28) 
(29) 
(30) 



We pre-compute the B matrices for a number of different incident mode frequencies and plasma-/? values. 
Finally, having constructed the transformation matrices T at a series of tube separations, we apply the 
standard MATLAB least-squares algorithm (backslash command) to solve equation ([50]) for Ai and then for 
A2. 



4. RESULTS 

First we consider the interaction of a pair of identical flux tubes at different separations with incident 
waves. The angle between the tubes is set to 712 = 0, meaning that the tubes lie on the x-axes of the 
coordinate systems of each other (see Figure [T]) . Undoubtedly there will be large changes in the scattering 
coefficients at different angles; for now we stick to this simple case, and leave further exploration of the 
parameter space to future endeavors. A pictorial representation of an incoming /-mode with respect to the 
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flux tubes 1,2 is shown in Figure [2] The resuhs of these interactions are shown in Figure [3l A number of 
conchisions may be drawn: (1) the strongest couphng is between the / mode and the flux tube; there is a 
rapid faU off with increasing radial order, (2) the length scale of the region occupied by the near-field region 
well approximated by 7r/fc^J, where 27r/fcP is the horizontal wave length of the incident p„ mode, (3) the 
scattering coefficients attain large values at small fiux tube separations; this has a number of consequences 
for wave absorption in plage, (4) the phases of the scattered waves exhibit highly unpredictable trends, even 
more so than the magnitudes of the scattering coefficient (Figure |4|), and (5) higher-/? flux tubes have more 
extended jackets than their lower- /3 counterparts; this may be attributed to the fact that when matching 
tube displacement eigenfunctions to the external waves, the jacket modes play a more significant role in the 
former case. The last point is further illustrated by the upper two panels of Figure [5] which show the isolated 
body displacements (,± of the tubes at plasma-/? = 0.1, 1. The fact that a larger number of kinks are seen 
in the higher-/? tube means a larger number of jacket modes must crowd the near field region, resulting in 
stronger coupling between the tubes. The lower two panels of Figure [S] demonstrate the stark differences 
between the isolated body and near-field coupled tube displacements; the structure, the number of nodes, 
and amplitude of the eigenfunction are seen to greatly increase when the jackets of the tubes are able to 
communicate. This is the very premise of multiple scattering. 

We also briefly investigate the interactional behavior of a pair of non-identical flux tubes in Figure [6l 
Because the magnetic fiux in each tube is held constant (see ^2.ip . the two tubes which have differing plasma- 
(3 of 0.1,1 also are of different radii, commensurate with equation ([8]). The scattering coefficient trends are 
unremarkable, showing differences in structure from those of Figure [3] but do not possess any noticeable 
features. Despite differences in the flux tube geometries, the coupling remains strong and continues to adhere 
to the half- wavelength rule-of-thumb near-field dimension. Lastly, in Figure [7l we graph the components of 
the wave field with the upper panels showing the tubes strongly interacting (a separation distance of 1.085 
Mm) while the isolated body case is displayed in the lower panel. In the strong interaction case, the near-field 
lobes of the two tubes are seen to be in communication, the scattering in one tube significantly influencing 
the other through this medium. The far-field in both cases are outward spirals, transporting some of the 
incident mode energy away from the scene of scattering. 



CONCLUSIONS 



The work of Bogdan fc ZweibeJ (1987) was among the first efforts to characterize and study the multiple 
scattering of waves by ffux tubes. Since then, f urther activity in this area was restri cted to the study of mul- 
tiple scattering by unstratified flux tubes (e.g. iBoedan fc Foxlll991t iKeppens et al.lll994) and su bseauentlv. 
single scattering by gravitationally stratified tubes (e.g. Bogdan et al.lll99?Hanasoge et ahlboosh . Including 
stratification when attempting to solve the full multiple scattering problem introduces manifold difficulties, 
especially witho ut the right set of tech niques . However, much progress in this regard and the availability of 
the technique of Kagemoto fc Yue (e.g. 19861 ) has allowed us to attempt this problem. 



There has been the general thinking that not only does gravity inhibit resonant absorption but may in 
fact prevent strong interactions between closely spaced flux tubes. The reason for this is the disparity in the 
depth structure between the eigenfunctions of external modes and the tube displacement; gravity introduces 
strong structural differences between these two quantities, almost certainly ruling out a strong and direct 
matching between the two. Thus resonant absorption in stratifled magnetized environments may be largely 
ruled out; but what about multiple scattering and spikes in scattering cross sections and phases? Is the 
theory of single scattering sufficient to describe the interaction of waves with clustered magnetic elements 
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Incident m=1,4mHz f mode 




X [Mm] 

Fig. 2. — Example of an incident m = 1 wave: ^'inc = iJi{kQr)e^^ . Tube 1 always sees the incident wave, 
whereas with increasing separation, tube 2 sees less and less of it. Notice that the presence of the second 
tube kills the m = ±1 symmetry when the tubes are close to each other. The symmetry is restored for wide 
tube separations, at which point, both tubes are practically isolated bodies. 



in the Sun? Certainly not, as demonstrated in fJD Fairly dramatic changes in the scattering coefficients are 
observed at close separations, of the order of several hundred kilometers, not unlike the distances between 
flux tubes in plage. The fact that scatter by a pair non-identical flux tubes also exhibits a similar trend 
(Figure [6|) demonstrates the robustness of multiple scattering type interactions. The scattering coefflcient is 
a proxy for the degree of absorption and mode mixing exhibited by the system; if nothing else, the jumps 
in the coefficients point to a loss of coherence of the incident modes, introducing a mechanism for wave 
damping. 



The coefficients behave in a quirky manner, rather similar to those obtained bv iBogdan fc Foxl (jl99ll ). 
Thus drawing out larger behavioral properties from these interactions is somewhat difficult without a more 
detailed search of the parameter space. Part of this quirkiness may be attributed to the interference between 
the sizable numbers of modes competing in the excitation of t he flux tubes. Through analyses of observations 
of thousands of small magnetic elements. iDuvall et al.l (j2006l ) succeeded in estimating the detailed scattering 
properties of these fe atures, concludi n g tha t the m = ±1 waves couple quite strongly with the magnetic 
tubes. Subsequently, iHanasoge et al.l (|2008r ) modeled these measurements in the single scattering limit in 



an attempt to constrain properties of the average magnetic element. However, these elements consist of a 
number of thin flux tubes, the scattering presumably in the strongly interacting regime due to their proximity. 
Merely with two tubes at a sequence of separations, the scattering coeflicients and phases display remarkably 
intricate behavior; extending this to a number of flux tubes at various random locations is arguably a difficult 
task. 

The near field in Figure [7] still remains to be directly detected in observations. One possible route 
towards this goal is to look for statistically significant auto-correlation signals in the vicinity of these small 
magnetic elements. This follows from the spatially stationary nature of the near-field modes; they merely 
pulse at the frequency of the incident wave in a thin envelop around the scatterer. However it is unclear 
how to use this information in a productive manner; evidently a greater understanding of these evanescent 
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Fig. 3. — Mode-mixing and scattering computed for / andpi incident waves at two diflFerent frequencies and 
azimuthal orders m = ±1 for a pair of flux tubes at a number of separations. These are the coefficients written 
according to the co-ordinate system centered around tube 1, the tube that encounters the untransformed 
zeroth order incident mode. Note that the symmetry of the m = ±1 scattering coefficients are destroyed 
due to the presence of tube 2; this symmetry is regained when the separation becomes large enough that 
tube 1 can be considered an isolated body. Clear signs of multiple scattering are observed for separations 
comparable to about half the horizontal wave length at the surface (~ 7r/fc^) of the incident waves, which are 
11.04 and 2.76 Mm for the / modes and 25.76 and 6.44 Mm for the pi modes at 2 and 4 mHz respectively. 
The degree of scatter exhibits a complex behavior at small separations and more importantly, shows rather 
large deviations from the isolated body values. Note also the differences between the m = ±1 modes; the 
presence of the second flux tube destroys any symmetries with respect to these two m waves. Also of interest 
is the fact that the /3 = 1 flux tubes couple via the near-field more strongly than the low-/? tubes; this is 
presumably because the displacement eigenfunction exhibits more depth-structure in the former, requiring 
a greater participation of the near- field modes, and thereby resulting in stronger coupling (see Figure [5]). 



waves is needed. 



The question relating to sunspot structure, as to whether one can expect to be able to helioseismically 
discern a monolith from a jelly fish, still shows promise, for it would appear that multiple scattering, despite 
the presence of strong gravitational stratific ation, still play s a m ajor role. Our preliminary conclusions, 
derived from Figures [3] and H] echo those of IKeppens et al.l (jl994r ) , who suggested that ensembles of flux 
tubes can absorb quite effectively, but not cause coherent phasing of the scattered waves. 



This work represents a small first step towards comprehending the action of multiple scattering in 
stratified magnetized environments. Much work remains to be done in terms of characterizing and under- 
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Fig. 4. — Phases of the scattered / mode. The parameter space is similar to that explored in Figure [S] No 
clear pattern is seen, indicating that interpreting the phases in a multiply scattered wave field is quite a 
non-trivial affair. 

standing the interactions of clusters of randomly located flux tubes and the impact on mode linewidths and 
observations of scattering in plage. 

The idea for this work occurred over conversations at the CSPA, Monash University, where S. M. H. was 
a visiting academic in 2008. Much of the computation was performed on the Stanford University solar group 
machines; thanks to Phil Scherrer for the use of these resources. Some part of this work was accomplished 
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THE T and B MATRICES 



The procedure outlined here is adapted from Kagemoto fc Yue ( 19861) . In order to construct T^;, we 
have the following formula for each |to, m'| < 1, n: 

[TaU = e^^— ')^-i/i^i,„,(fc„i?.0 (Al) 
for n < np, where c = Sn + m' + 2 and d — 3n + m + 2, and 

{Tu)cd = e'(™-'"')^"(-l)'"'''^m-m'(fcni?.i) (A2) 

for n > np, where c,d are as above. The means of constructing Bj are as follows. For every tuplet 
(|m| < l,n) an incident wave with the appropriate eigenfunction and radial behaviour of unit amplitude is 
chosen. If ri < rip, the incident mode is propagating and has a Jm(fcPr) type horizontal behavior, whereas if 
n > np, the horizontal part of the eigenfunction is given by /^(fc^r). For each incident mode, the scattering 
into all other modes (at constant frequency and m) is computed. Having thus computed the amplitudes of 
all the scattered waves for each incident mode, we fill up the matrix Bj as follows: 

(B,)cd = (A3) 

where /3 is the scatter into the mode (m, n') by incident wave (m, n). The indices are c — 3n' + m + 2, d = 
3ri + m + 2. We only consider |r7i| < 1 since the tube is insensitive to other azimuthal orders in this theory. 
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Normalized tube displacement eigenfunctions for a 4 mHz, / incident mode 
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Fig. 5. — Real parts of the normalized tube displacement eigenfunction (Ci/fcfJ for flux tube 1 for /3 = 0.1, 1, 
and separations, A = 1.20,3.25 Mm. At this frequency, tubes at a separation distance of 3.25 Mm can be 
considered isolated bodies. The first aspect take note of is that the displacement eigenfunction of the 
isolated body /3 = 0.1 tube has fewer nodes and therefore less structure in depth than the (3 — 1 tube (upper 
two panels). Consequently, a greater near- field participation is required to successfully match the [3—1 
tube displacements with the external wave eigenfunctions, resulting in a more diverse near-field than the 
P — Q.l case. Secondly, the lower two panels show how complicated the tube displacement eigenfunction 
becomes when the near-fields of the two tubes are in communication. Also, it is important to note that the 
displacement amplitudes in the strong interaction case are larger by a factor of 10 or more than the isolated 
body counterpart. This aspect underscores the very premise of multiple scattering. 



B. TRANSFORMATION OF THE INCIDENT WAVE FIELD 



We use Graf's addition formula fe.g. lAbramowitz fc Stegunlll964[) : 

oc 



(Bl) 
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p -f scatter, v = 4 mHz 
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Fig. 6. — Scattering coefficients of the (3—1 tube for an interacting pair of non-identical (3 = 0.1, 1 tubes. 
The behavior is unremarkable, not greatly differing from that in Figure [S] Note that because the magnetic 
flux in the tube is held constant, changes in [3 necessarily imply a change in flux tube radius. The incident 
/ and pi modes have horizontal wavelengths of 2.76 and 6.44 Mm respectively. 



An m — Q wave seen at the axis center of tube i produces the following components at flux tube I whose 
center is located at distance Ra (the m = [—1, 0, 1] waves): 



( Ji(fc„i?,i)e'^" MknRa) J-i(fc„i?,i)e^'^")- 
Similarly m = ±1 waves seen at tube i produce at I: 

[MknRu) J-i{k„Ru)e-'''^' J-2{knRa)e'^'^^' ) , 



(B2) 



(B3) 
(B4) 



respectively. Since in our theory, the tube is insensitive to |m| > 1, we do not include these terms in the 
transformation. 
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